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Abstract 



We present a new program package for calculating one-loop Feynman inte- 
grals, based on a new method avoiding Feynman parametrization and the con- 
traction due to Passarino and Veltman. The package is calculating one-, two- 
and three-point functions both algebraically and numerically to all tensor cases. 
This program is written as a package for Maple. An additional Mathematica 
version is planned later. 
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PROGRAM SUMMARY 

Title of program: ONELOOP.MA 

Computer: (i) VAX 4000-90 (ii) PC 486 DX 

Operating system: (i) VMS (ii) MS-DOS 6.0 

Programming language: MAPLE V 

Subprograms used: CFCN.MA, R.MA, SIMPLE.MA, PV.MA 

No. of lines in distributed program, including subprograms: 1650 

Keywords: Electroweak theory, Feynman diagrams, one-loop corrections, renormalization 

Nature of physical problem: The theoretical determination of cross sections in particle pro- 
cesses requires the calculation of radiative corrections. The most important contribution 
comes from the level of one-loop Feynman diagrams which arise from the model under con- 
sideration, usually the standard model of elementary particles 

Method of solution: This package is designed to evaluate automatically one-loop integrals. 
It is making use of the properties of TZ functions, a class of special functions which simplifies 
the evaluation of Feynman integrals 

Typical running time: The algebraic evaluation of a two-point function up to the rank 2 
tensor on a VAX 4000-90 needs less than 20s 

LONG WRITE-UP 
1 Introduction 

In particle physics the calculation of one- and higher loop corrections to particle pro- 
cesses is necessary to keep track with the increasing accuracy of particle colliders. In 
this context it is very useful if a fast and automatized calculation of the contributing 
Feynman diagramms is provided by a computer. On the most basic loop correction 
order, the one-loop level, several attempts for such programs are already available, 



like FF |13| - calculating scalar one- loop integrals in Fortran - or FeynCalc |12] - a 



Mathematica |11] package - which exhibits the usual Passarino-Veltman coefficients 



of tensor integrals in terms of scalar integrals. Both programs are using the stan- 
dard procedure for calculating one-loop diagrams, namely Feynman parametrization, 
Wick rotation and the Passarino-Veltman contraction for tensor structure (cf. [^ and 

§)■ 

Although this standard procedure became very popular during the last decade, 
it turned out to be rather inconvenient for being implemented on a computer if the 
integral involves higher tensor structure. For that reason we investigated in a new 
method for calculating arbitrary one-loop integrals which is already published for 
the two- and three-point case (cf. Q and [Q]). 

Our approach is completely different from the standard procedure and the struc- 



ture of other packages. We are not forced to invert big matrices which is the most 
difficult step in computing the tensor structure of the integrals. In our method the 
loop momentum integrations are performed directly, using the residue theorem and 
introducing TZ functions (cf. [^), a class of special functions with some nice and use- 
ful features. Instead of big matrices we get some recursion relations for the involved 
TZ functions (cf. Q or |p ) , which allow a fast and memory efficient evaluation of the 
integrals under consideration. 

The program is designed to reduce the TZ functions to a set of fundamental 
functions which are rewritten in terms of logarithms and dilogarithms. Up to now 
the different one-, two- and three-point functions are programmed as procedures in 
Maple [10|. Therefore our results are available numerically as well as algebraically. 
We tested our procedures using Maple V on a VAX 4000-90 workstation. 



2 Description of the method 

A detailed description of the method may be found in Q or Q. Here we would 
like to give a brief introductory overview only. The calculation is based on the well- 
known dimensional regularization method, as described for example in g. Instead of 
introducing Feynman parameters the integral is split up into orthogonal and parallel 
space (cf. [P), so that for example the two-point integral takes on the form 



5(P"Pi)(g,mi,m2) 



(1) 



D-l 



dh 



dl^ if-^ 



ihy'ih) 



pi 



[(^11 + q)"^ -l\-m\ + ig] [E -l\-ml + ig] 



Here we present the most general case, depending on one external momentum q 
and two arbitrary masses mi and m2- The splitting of space requires a distinction 
between l\\ and l±_ which we define as 

l-Q 



u 



li-f 



(2) 



In this notation l\\ describes the component of / which is parallel to the external 
momentum g, whereas l±_ represents the orthogonal complement. 

The integrals (||) can be reduced to a sum of TZ functions by using ( A.l ) with the 
help of partial fraction and residue theorem. The complete procedure is described 
in detail in and |]l|. The result in the two-point-case may be written according to 
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where e = ^2~ represents the usual dimensional regularization parameter. We used 
the abbreviations 

2 2 

mi — mi 

A = —^ 

2q 

r , _ J a, a even 

^""^ ~ 1 a + 1, a odd 

The result of the three-point integral is obtained by a quite similar procedure and 
may be found in Q or [Q]. 

3 Evaluation and expansion of the IZ functions 

After having expressed all possible integrals in terms of TZ functions we get easy and 
fast computable results by applying some recurrence relations for the TZ functions 
under consideration. 

We achieve these relations by first reducing the number of different TZ functions 
to a set of fundamental TZ functions. This is possible due to the already mentioned 
features of the TZ functions, especially the relation ( [A.2| ) and ( |A.3| ). The reduction 
takes place in two steps. First of all we increase the parameter of the TZ function 
to get it close to or 1. This can easily be achieved by using formula ( [A.2| ). In the 
second step we reduce the index of the TZ function by virtue of a convenient formula 
derived from ( |A.3| ). In the case of the two point-function it has the form: 

1Zt{bi,b2;zi,Z2) = ., n_. [{{t- 1)^1^2) TZt^2{biM\ zi,Z2) (3) 

+ {{l-t- bi)zi + {l-t- 62)^2) TZt-i{biM; ^1, ^2)] 
Now the power of this method is obvious. Both steps mentioned above may be 



performed easily by computer, due to the recursive definition of the formulae (A. 2) 
and O. 

Up to now we have acquired a result which is exact to all powers of e, though 
we are interested only in an expansion in powers of e. Therefore we now substitute 
the expansion for the fundamental TZ functions in terms of e, resulting in simply 
computable functions like logarithms and dilogarithms. This is a rather cumbersome 
task but has to be performed only once for each N-point function. The procedure, 
especially its subtleties concerning analytical continuation, is described in detail in 
another paper Q. For completeness we here just cite the expansion for the scalar 
two-point case 
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The three-point integrals 
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are calculated in a quite similar way (cf. Q or [^). Again we treat the general case 
with two independent external momenta qi,q2 and three masses mi, ?7i2, "23. For the 
momentum components we write in a quite natural generalization: 
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The evaluation is resulting in the fundamental 7^ functions 
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noting that our relations are valid in all kinematical regions, so that the results of 
apply to all the cases of our program. For a detailed discussion see Q. 



4 The one-loop functions 

Our aim is to calculate one-loop integrals both algebraically and numerically. There- 
fore the algorithm is implemented in Maple jl^, a 'language' for symbolic compu- 
tations. To make the implementation for any user most comfortable, the procedures 



are bound to a package, which may be read from the Maple system during run-time. 
The user-interface consists of three elementary functions 

• OneLooplPt(p, m) 

• OneLoop2Pt{po,pi,q,mi,m2) 

• OneLoop3Pt(po,Pi,P2, ^i, 920, ^21, "^i, "^2, "^s) 

The functions directly correspond to the following integrals 

• One-point function: 



A(p){ 



m] 



d'^l 



{ly 



\p — m? + ig] 



Two-point function: 



D-l 



n 



dill Idli l?-^ 



(^iirG±) 



pi 



[(^11 + Qf -11-1711 + ig] [Ij -l\-ml+ ig] 



Three-point function: 

C(PoPiP2)(g^^^2,mi,m2,m3) 



D-2 OO 

27r^- 



OO 



r(^) 



dL 



oil / "HI 

OO — OO 



dUn Idli l?-^ 



WHhwyHh)'' 



[{low + ^i)^ - ^?|| -ll-mj + ig] 



[{low + «20)^ - (^l|| + 921)2 _ /^ _ ^2 ^ ig^ [/2 ^ _ /2 ^ _l2^_ml+ ig] 

Again we follow our notation which distinguishes between parallel and orthogonal 
space. The masses are denoted by rrii and the definitions of (|2|) and (^) for the 
momenta are of course still valid. We would like to emphasize the fact that in our 
notation the indices pi represent the powers of the different components of the loop 
momentum I and should not be mixed with Lorentz indices of the corresponding 
tensor. This notation is perhaps not very often used, but it corresponds directly to 
our method of integration. The usual notation may be recovered by an additional 
algorithm described in section |5|. 

The results of the one-loop functions are Laurent expansions in terms of the ul- 
traviolet regulator e. The output of the described procedures consists of a list where 
the significant coefficients of this expansion are denoted. 



input 


output 


OneLooplPt(p, m) 


[e^^-term, e^'-term] 


OneLoop2Pt(pO)Pi) 9, mi, 1112) 


[£~^-term, e^'-term] 


OneLoop3Pt {po,Pi,P2,qi, 

q20,q2i,rni, 7712,1713) 


[£~'^-term, £~^-term, 

£^-term, {abb.}] 



In the output "e~^-terin" is meant to represent the divergent part (the coefficient 
of -) whereas "e^-term" is describing the finite part. Of course the e~^-terin has 
to vanish in any case, but is left in order to check the accuracy of the program in 
numerical calculations. In the three-point case "abb." is written for some abbrevi- 
ations which are introduced to get a more compact result. For the same reason we 
define a function A which collects the ocurring rj functions in the three-point case: 



A{x,y) =ln{y) 
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(6) 



Up to now the program is not regulating the infrared divergence of the three-point 
function - which may occur if one mass is set to - by virtue of an infrared dimension 
parameter eir. Therefore the user is expected to avoid this kind of divergence by 
introducing a (small) mass parameter instead of 0. 

Another difficulty arises if the external momenta are collinear, which means that 
521 = 0. Although this case is simpler to calculate, because the parallel space is only 
one-dimensional, the program is not taught to deal with this special case up to now. 
The user is asked to use a small g2i-component instead. 

As long as the arguments of the OneLoop functions are symbols the output is 
given algebraically, whereas of course numbers inserted for the arguments imply a 
numerical result. Moreover, in the numerical case, the program sets the value of g, 
the imaginary part of the propagators, five digits higher than the numerical accuracy 
of the whole calculation. In the algebraic case just 'rho' is returned. 

For different calculations of the same OneLoop function it is of course rather 
inconvenient to start the whole program again. Therefore our package includes 
the functions OneLoopLiblPt(n), OneLoopLib2Pt(n) and OneLoopLib3Pt(n) which 
generate a library of all OneLooplPt, OneLoop2Pt and OneLoopSPt functions re- 
spectively up to the tensor rank n. Each function is written to one ffie. The ffienames 
correspond to the integrals in the following way: 



ONEO.RES 

TWOOO.RES 

TWOIO.RES 

THREEOOO.RES 

THREEIOO.RES 

THREEOIO.RES 



OneLooplPt(0, m) 
OneLoop2Pt(0, 0, q, mi , 7712) 
OneLoop2Pt(l, 0, q, mi , 7712) 
OneLoop3Pt(0, 0, 0, gi, 920, 921, w-i, "1-2, "t-s) 
OneLoop3Pt(l, 0, 0, gi, 920, 921, w-i, ■m2, ms) 
OneLoop3Pt(0, 1, 0, gi, 920, 921, "t-i, "t-2, "t-s) 



To specify the path of the library the variable "tensorpath" may be assigned to the 
desired directory. 

5 The Passarino-Veltman coefficients 



To receive a more familiar representation of the output a procedure is added 
which extracts the Passarino-Veltmann coefficients from our programs, PassVeltlPt, 
PassVelt2Pt and PassVeltSPt for the different n-point cases. The following types of 
arguments - here demonstrated for instance for the two-point function - are possible: 



• PassVelt2Pt(n): The tensor decomposition of the rank M tensor two-point 
function 

[{I + qY — m\ + ig] [P —7712+ ^^] 

is given. The procedure returns a hst consisting in the different coefficients, 
which are expressed in terms of the OneLoop tensor integrals, and the defining 
equation for the coefficients, for instance in the case M = 2: 

■ OneLoop2Pt(0, 2) OneLoop2Pt(0, 2) 

3 - 2e q^ (—3 + 2e) 

OneLoop2Pt(2, 0) 

-I 2 ' ^209mi9m2 + ^215^1^2 

• PassVelt2Pt(n,full): The function returns the same as above, but inserts the 
results of the OneLoop tensor integrals explicitly. 

• PassVelt2Pt(n, g, mi,r?T,2): The 
function returns the same as PassVelt2Pt(n,full), but expressed in the user 
defined terms for q,mi and m2, using not the default parameters. Of course 
numerical values are also possible. 

In the three-point case a similar decomposition of the rank M tensor three-point 
function 



[{I + qi)'^ - mf + ig] [(/ + ^2)^ -ml + ig] [P -ml + ig] 

is returned 

If the variable "tensorpath" is set to some path, the PassVelt procedure assumes 
that the OneLoop functions are written there as described in the previous section. 
It then reads the necessary integrals from this directory. Otherwise, if "tensorpath" 
is not assigned, PassVelt is forced to calculate the necessary OneLoop functions by 
its own. 



6 An example 

To demonstrate how our package works we illustrate the calculation of the electron 
self-energy in QED. Starting with QED Feynman rules we get: 
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Now we split into orthogonal and parallel space noting that the 7 matrix vector is 
split into 7|| and 7_l: 
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At this stage of calculation we now take the opportunity to incorporate our 
procedures which solve the ocurring integrals iJ^W'Pi-'. To be specific we have 
to compute the functions OneLoop2Pt(l, 0, g, m, 0), OneLoop2Pt(0, 1, g, m, 0) and 
OneLoop2Pt(0,0,g,m,0): 
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Inserting these functions and simplifying the expression one may rewrite the result 
in the following way: 
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which is the well-known result for the electron self-energy. 



+ 0{e) 



(8) 



7 Conclusion and further projects 

Our aim was to introduce a program package that gives any computer user the pos- 
sibilty to calculate one-loop particle processes without much effort in manpower and 
hardware. At this stage of development the program is understood to be the first 
step of a larger project. We do not claim that the package is up to now completely 
bug-free and optimized. In fact we are planning to incorporate the possibility to reg- 
ulate infrared divergences also in the dimensional regularization scheme and taking 
into account higher powers of propagators as well as the special case of collinear 
external momenta. 

As a further project we plan to extend this package to the four- and five-point 
function. Moreover we want to export our program package to Mathematica |11], 
because from our point of view Mathematica at the moment is the more common 
language than Maple. Planned is also a menu-driven user interface. But this might 
be rather difficult, because we haven't yet found a convenient C-language package 
for algebraic calculations. (For hints contact L. Briicher) 
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Appendix 



Conventions 

D = 4: — 2e dimension in the dimensional regularisation method 
TZt{b;z) TZ function, t is cahed index, 

h = (61, ... , hn) parameters, z = (zi, . . . , z„) arguments 
(3 = Yll=i bi sum of parameters of an 7^ function 

Some fundamental properties of TZ functions 

• formula for integrals 

/x"~^ flizi + Wix)-^' dx = B{P - a,a)na^p{h, —) f[w^''' (A.l) 

• increasing of parameters 

/37^i(b,z) = (/? + i)7^i(b + ei,z) - izi7^i_l(b + ei,z) (A.2) 

• associated functions (k is the number of parameters) 

k 

J2 A\''\t,h,z)n-t-iih,z) = (A.3) 

i=0 

with 

E^ (z) is the elementary symmetric function, gained from 

k 

il+xzi)---{l + xzk)=J2 2^'^f ^ (-2) (A-5) 

1=0 
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